Competition of Zener and polaron phases in doped CMR manganites 
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Inspired by the strong experimental evidence for the coexistence of localized and itinerant charge 
carriers close to the metal-insulator transition in the ferromagnetic phase of colossal magnetore- 
sistive manganese perovskites, for a theoretical description of the CMR transition we propose a 
two-phase scenario with percolative characteristics between equal-density polaron and Zener band- 
electron states. We find that the subtle balance between these two states with distinctly different 
electronic properties can be readily influenced by varying physical parameters, producing various 
"colossal" effects, such as the large magnetization and conductivity changes in the vicinity of the 
transition temperature. 

PACS numbers: 71.10.-w, 75.30.Vn, 71.38.Ht, 75.30.Kz 
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I. INTRODUCTION 

The transition from a metallic ferromagnetic low- 
temperature phase to an insulating paramagnetic 
high-temperature phase observed in some hole-doped 
manganese oxides (such as the perovskite family 
Lai_ x [Sr, Ca] x Mn03) is associated with an unusual dra- 
matic change in their electronic and magnetic properties. 
This includes a spectacularly large negative magnetore- 
sistive response to an applied magnetic field - sometimes 
termed colossal magnetoresistance (CMR) - which might 
have important technological applications [Q. Starting 
with the pioneer papers of Jonker and van Santen half a 
century ago 0, this challenging behavior has stimulated 
a considerable amount of experimental and theoretical 
work |^-^|, however, even much of the basic physics of 
the CMR still remains controversial. 

Early studies on lanthanum manganites concentrated 
on the link between magnetic correlations and transport, 
and attributed the low-T metallic behavior to Zener's 
double-exchange mechanism KH§], which maximizes the 



hopping of a strongly Hund's rule coupled Mn e g -electron 
in a polarized background of the core spins (Mn ti g - 
electrons). The quantum version of this process has been 
described by Kubo and Ohata ||. Although there is 
no much controversy about the qualitative validity of 
the double-exchange scenario to stabilize a ferromagnetic 
state, it has been argued that physics beyond double- 
exchange is important not only to explain the very com- 
plex phase diagram of the m ang anites pO| , pl| but also 
the CMR transition itself HQ. The difficulty is that 
magnetic scattering of itinerant charge carriers (doped 
holes) from enhanced fluctuations near T c is the exclu- 
sive mechanism to drive the metal-insulator transition in 
the double-exchange-only models. However, even com- 
plete spin disorder does not lead to a significant reduc- 
tion of the electronic bandwidth, and therefore cannot 
account for the observed scattering rate p4]dT6|] (cf. also 
the discussion in Refs. ]j~7| , |l8|] ). 



In view of this problem, it has been suggested that or- 
bital 19-23[j and, in particular, lattice effects [fl2 24-2^1 
are crucial in explaining the CMR phenomenon. The ar- 
gument was that the electron-lattice coupling is known 
to be strong in at least some members of the perovskite 
manganese family and tends to localize doped hole carri- 
ers as small polarons by changing the ratio of the energy 
gain through polaron formation to the conduction band 
kinetic energy 14,2l|]. Clearly this ratio is extremely sen- 
sitive to changes of the magnetic correlations by vary- 
ing magnetic field, temperature and carrier concentra- 
tion |lj,|2^]. There are two types of lattice distortions 
which are important in manganites. First the partially 
filled e g states of the Mn 3+ ion are Jahn- Teller active, 
i.e., the system can gain energy from a quadrupolar sym- 
metry elongation of the oxygen octahedra which lifts the 



e g degeneracy 



A second possible deformation 



is an isotropic shrinking of a MnOg octahedron. This 
"breathing" -type distortion couples to changes in the e g 
charge density, i.e., is always associated with the presence 
of an Mn 4+ ion. In the lightly doped region of the phase 
diagram, holes are pinned onto Mn 4+ sites by locally "un- 
doing" the cooperative Jahn- Teller distortion, i.e., form- 
ing so-called "anti Jahn- Teller polarons" [Bifl. As pointed 
out in a recent paper by Billinge et al. |32|], this per- 
spective is not appropriate, or at least misleading, in the 
heavily doped material, when the high-temperature pola- 
ronic state is approached from the ferromagnetic metal- 
lic (Zener) state, because there are initially almost no 
Jahn- Teller distorted octahedra. In this regime, both 
breathing- mode collapsed (Mn 4+ ) and Jahn- Teller dis- 
torted (Mn 3+ ) sites are created simultaneously when the 
holes are localized in passing the metal-insulator transi- 
tion. To avoid any confusion, in what follows, we refer 
to a "polaron" as a doped charge carrier (hole) which is 
quasi-localized with an associated lattice distortion |}2j . 
The relevance of small polaron transport above T c is obvi- 
ous from the activated behavior of the conductivity |33| . 
Consequently many theoretical studies focused on po- 
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laronic approaches |1^J3J |39j] . Polaronic features have 
been established by a variety of experiments. For exam- 
ple, high-temperature thermopower jl(],fy] and Hall mo- 
bility measurements [ |42| confirmed the polaronic nature 
of charge carriers in the paramagnetic phase. More di- 
rectly the existence of polarons has been demonstrated by 
atomic pair distribution p3[ , X-ray and neutron scatter- 
ing studies . Interestingly it seems that the charge 
carriers partly retain their polaronic character well be- 
low T c , as proved, e.g., by neutron pair-distribution- 
function analysis p7| and very recent resistivity measure- 
ments |4^]. Moreover, there have been predictions, based 
on XAFS data |^9| , that small octahedral distortions per- 
sist at low temperature, forming a nonuniform metallic 
state. Particularly striking in this respect is the direct 
relationship between the structural distortions and the 
magnetism of the CMR perovskites |30|. That means, 
even the nature of the ferromagnetic low-temperature 
phase is not yet completely resolved. Obviously man- 
ganese oxides, above and below T c , are in the subtle 
regime where many different tendencies are in compe- 
tition and it seems that more refined ideas are needed to 
explain the main properties of these materials. 

Realizing that intrinsic inhomogeneities and mixed- 
phase characteristics exist and might play a key role 
in manganites, two-phase models describing the coex- 
istence and interplay of localized and itinerant carriers 
have recently attracted a lot of attention (for a review 
see Ref. ||). For example, phase separation scenarios 
involving phases with different electronic densities have 
been adopted to describe the mixed-phase tendencies in 
manganites |^3|,|lJ. These approaches are particularly 
meaningful in the limits of small (a; < 0.1) and high 
{x ~ 1) hole densities, where nanometer scale coexist- 
ing clusters have been reported. In the CMR regime 
(0.15 < x < 0.5), however, several experiments re- 
ported clusters as large as micrometers in size |32],[33| . Of 
course, such /mi-sized domains, if charged, are energeti- 
cally unstable because of the electro-neutrality condition 
enforced by long-range Coulomb repulsion. Therefore 
an alternative concept, where the metal insulator tran- 
sition and the associated magnetoresistance behaviour is 
viewed as a percolation phenomenon, was analyzed theo- 
retically [p^-|56| , and now has a substantial experimental 
support by pulsed neutron J57| , electron microscopy fl53[ , 
Mossbauer |5f| and scanning tunnel spectroscopy ]5^ ] 
studies. 

Motivated by this situation, in the present work we 
propose a novel approach to the metal-insulator tran- 
sition in manganites which takes into account the per- 
colative coexistence of two "intertwined" equal-density 
phases: metallic double-exchange dominated and pola- 
ronic insulating. The transition is driven by a feedback 
effect which, at T Cl abruptly lowers the number of de- 
localized holes, i.e., conducting sites, leading to an col- 
lapse of the bandwidth of the Zener state. The physical 
picture behind our approach is corroborated by X-ray- 
absorption fine structure M and pair distribution |32] 



data, which indicate that charge localized and delocalized 
phases coexist close to the CMR transition ]50],[52| . Fur- 
ther support comes from zero-field muon spin relaxation 
and neutron spin echo measurements |59| l , reporting two 
time scales in the ferromagnetic phase of Lai_ x Sr x Mn03, 
where the more rapidly relaxing component has been at- 
tributed to spins inside the overlapping metallic region 
and the slower Mn-ion relaxation rate was associated 
with polarons which occupy a diminishing volume frac- 
tion as the temperature is lowered. 



II. TWO-PHASE MODEL FOR THE CMR 
TRANSITION 

To set up a two-component description of the metal- 
insulator transition in CMR manganites, let us first char- 
acterize the two competing phases. 



A. Delocalized Zener state 

The degeneracy of e g states in (cubic) manganese ox- 
ides implies that the itinerant e g charge carriers carry an 
orbital degree of freedom (put in mind that the Jahn- 
Teller distortion lifting the e g degeneracy vanishes in the 
highly-doped low-temperature metallic phase). Conse- 
quently the e g electron transfer amplitudes between two 
Mn sites depends on the orbital orientation. Within 
an {9 = |3z 2 — r 2 ), e = \x 2 — y 2 )} orbital basis, 
the anisotropic hopping matrix elements between bonds 
along the x, y, and z directions are given by pO||jl[ 
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(a, (5 € {9, e} denote orbital indices). We note that the 
orbital pseudospin is not a conserved quantity. The dis- 
persion of the corresponding (noninteracting) two-band 
(C = ±) tight-binding Hamiltonian is 
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The band structure (g) yields the density of states (DOS) 
depicted in Fig. |l] (cf. also Ref. ^8|). Although we em- 
ploy this DOS for all the calculations presented below, we 
would like to emphasize that the results will not change 
dramatically using an (isotropic) simple cubic or even 
constant DOS. 
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FIG. 1. Density of states corresponding to the band disper- 
sion (^) [W = 12t denotes the bandwidth]. Dotted lines give 
the DOS of the two ("-subbands. At optimal doping x ~ 0.3, 
the chemical potential = fJ-{T = 0) is located near the 
lower maximum of the DOS. 

Starting from an Kondo lattice Hamiltonian, in the 
limits of strong Hund's rule coupling and large on-site 
Hubbard interaction (U ^ Jh 3> which reflects 

the situation in the manganites), an effective transport 
Hamiltonian for (spinless) e g holes can be derived ||,|l6| . 
Note that because this Hamiltonian acts in a projected 
Hilbert space without double occupied sites, in Fig. |l| 
the integrated DOS is normalized to one, leading to an 
additional prefactor 1/2 in the (k, (^-summations. Then, 
combining the action of Zener's double-exchange with the 
percolative approach, the following effective hole band 
results 



(3) 



That is, in our model the renormalization of the band 
energy of the Zener state is driven by two mechanisms. 

At first, introducing an effective field A = (3g^BH^, s 
that tends to order the ion spins in z-direction, we have 
the temperature- and field-dependent band narrowing 
due to the Kubo-Ohata factor 
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Next, the percolative aspects of the metal-insulator 
transition discussed above imply that, at least just above 
x c (critical concentration for the occurence of the ferro- 
magnetic metallic state at T = 0) or below T c , there 
exist insulating (polaronic) enclaves sparsely embedded 
in the conducting ferromagnetic (Zener) phase. Thus the 
ferromagnetic phase occupies a volume smaller than the 
sample volume, i.e. < N, where N is the total num- 

ber of sites and denotes the number of ions in the 
FM phase. Of course the sizes and shapes of the insu- 
lating microscopic islands are difficult to access quanti- 
tatively. For simplicity, within our effective two-fluid ap- 
proach, we assume, that the hole hopping amplitude has 
the value t inside the conducting region (corresponding 
to the effective transport Hamiltonian (28) of previous 
work |l6|]) and zero elsewhere. Using Fourier analysis, 
the homogeneous component of this spatially dependent 
hopping amplitude, t — N N t, gives the renormalization 
of the bandwidth, whereas the long-wavelength Fourier- 
components cause the scattering of the carriers. Natu- 
rally the size of the ferromagnetic region, or equivalently 
the "probability" 
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has to be determined self-consistently (see Sec. II C). 
This introduces a feedback effect. Additional justifica- 
tion for the assumed reduction of the bandwidth propor- 
tional to the fraction of the ferromagnetic region can be 
obtained from the numerical calculation of the density of 
states outlined in the appendix for a simple tight-binding 
site percolation model without feedback. 



B. Localized polaronic state 

Let us denote the number of ions and localized elec- 
tron vacancies (holes) in the polaronic phase by 
and respectively. Of course, At(p) = N - Nut and 

N^ p) = Nh- hold, where N h is the total number 
of holes, i.e., x — N^/N . Then the energy gain due to 
the Jahn- Teller splitting on localized electron sites with- 
out the influence of vacancies is weakened according to 



(jy(p) _ N ( t j ?) )E 1 = {x- 1 - l)ExN^', where E x describes 
an effective Jahn- Teller energy in the polaronic regions. 
At the same time, if a doped hole is localized at a cer- 
tain site, forming a small polaron, a breathing distortion 



t(p) 



Here 5 = 5 + 1/2 refers to the total 3d spin of an Mn 3 + may occ ' r whkh lowerg the energy of the unoccupied 
ion when an e g electron is present, and S = 3/2 is the 
total core spin formed by the three <2g electrons. The 



factor M), yielding an effective hole transfer amplitude 



* = 7g [SA]t, 



(5) 



was obtained averaging the transfer matrix element for 
a single nearest-neighbour Mn 3+ -Mn 4+ bond over all 
values and directions of the total bond spin St (cf. 
Refs. |B). 



level by the familiar polaron shift, E p — —g mo (see, e.g., 
Rcf. |62j), relative to its energy in an ideal structure. 
In addition, a reduction of the Jahn- Teller distortions in 
the neighborhood of the localized electron vacancy takes 
place. Both effects add up to an effective polaron binding 
energy E2 < per hole. Neglecting the exponential small 
polaronic bandwidth, the polaronic phase, realized only 
m a fraction = jy(p) /jy of the sample, can be repre- 
sented approximately by spinless fermions (holes) having 
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the following site-independent energy 
£ P = (a^ 1 - 1) Ei + E 2 



(7) 



As follows from a fitting of the resistivity data to a small 
polaron model |53), in the region 0.2 < x < 0.5, the 
x-dependence of the parameter Ei is rather weak, and 
therefore has been omitted in (Q). 



C. Self-consistency equations 

The basic assumption of our mixed-phase scenario is 
the coexistence of metallic and insulating clusters with 
equal hole den sity, i.e ., in accordance with recent experi- 
mental results |5!^,|32j (see also Ref . 56 ) , we assume that 
there is no large-scale separation of Mn 3+ and Mn 4+ ions 
in the CMR doping regime. This means 
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but, of course, in general we have ^ and pW ^ 



Denoting by 1(H) the energy of the whole two-phase 
system within mean-field approximation and by S the 
corresponding entropy, T = (H.) — TS is an upper bound 
for the free energy. Introducing the grand-canonical po- 
tentials 



and 



n(/) = _ J_ y ln[l + e^-^) 
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(10) 



for holes in the ferromagnetic and polaronic phases, re- 
spectively, 



results, where 



(11) 
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represents the mean-field ion-spin entropy, and 

Vg[z] = sinh(z) coth(^) + cosh(z) , 
^W = % i coth(M±iz)- 5 Lcoth(^). 



(12) 

(13) 
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If we use, instead of the model (28), the spin-dependent 
transport Hamiltonian (29) of our previous work f pTj| , 
assuming that the correlations of the spin background 



change on a time scale large compared with the hole 
hopping frequency, i.e., we replace the effective hopping 
amplitude (0) by 



[S(1 + B S [SX})Y 
(2S)(25 + 1) 

the ion-spin entropy takes the form 
5 W =-.V 



t, 



(15) 



{p (f \lnu § [SX}-XSBs[SX}) + p (p) lai/ 5 [0]}. 

(16) 



Most significantly, at given temperature T and doping 
level x, both the ordering field for the ferromagnetic state 
(A) and the size of the Zener phase (N^\ or alternatively 

the hole fraction N^) have to be determined in a self- 
consistent way, minimizing the free energy ( |ll] ) on the 
hyperplane p,(X,N^) given by 
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After all, the magnetization (more exact the averaged 
spin z-component) can be calculated from 



M = (1 - x)S \p( f) B s [SX] + p ip) B s [0] 
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III. NUMERICAL RESULTS AND DISCUSSION 

The changes in magnetization as a function of tem- 
perature are shown in Fig. |^ for different doping lev- 
els. The efficiency of the feedback effect is quite obvi- 
ous. Omitting the p^-factor in Eq. (||), a rather smooth 
variation of M(T) results and the critical temperatures 
for the disappearance of the Zener state, T c (x), are too 
high as compared with the experimentally observed ones. 
Combining percolative and double-exchange mechanisms 
within our simple two-phase model, the magnetization 
exhibits a first-order transition, whereby T c is reduced 
substantially. It can be expected that the inclusion of 
a small but finite polaron bandwidth (instead of the lo- 
calized level (0)), softens the abrupt transition to some 
extent but in any case we will find a rather sharp drop 
of M (T) near T c . 

Naturally the bandwidth of the Zener state, displayed 
in Fig. H for a doping level x — 0.3, reflects the behavior 
of the magnetization. Without feedback we found a mod- 
erate band narrowing, whereas a radical shrinking occurs 
in the feedback model, which can be traced back to the 
collapse of the ferromagnetic metallic region at T c . The 
use of the spin-dependent hopping amplitude (|l5| ) leads 
to a modification of the Zener bandwidth (and /z) only 
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for the scenario without feedback. Here, temperature- 
independent band edges correspond to complete spin dis- 
order. Including thep^-factor, both models ( |l5| ) and (||) 
yield nearly identical results, which again indicates the 
dominance of the percolative feedback mechanism. 




400 600 
T[K] 

FIG. 2. Magnetization M, normalized by Mo = S — x/2, 
as a function of temperature T at various doping levels 
x — 0.175, . . . ,0.4. Results are shown for the models with 
(bold lines) and without (thin lines) feedback using parame- 
ters Ex = -0.125 eV and E 2 = -0.25 eV, and t = 0.3 eV. 
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FIG. 3. Temperature-dependences of the Zener band 
(shaded region) and of the positions of the polaronic level 
(e p ) and chemical potential (/i) without (a) and with (b) feed- 
back. Results are presented at x = 0.3 (Ei = —0.125 eV, 
E2 = —0.25 eV; t — 0.3 eV). Dashed lines display the band 
edges resulting from ( |l5| ) instead of (^). 
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FIG. 4. Phase diagram of the mixed-phase Zener-polaron 
model with feedback. The inset shows the fraction of the 
Zener phase as a function of temperature. 

Figure 4 presents the central result of our work, the 
x-T phase diagram of the percolative two-phase model. 
First of all we would like to state that it is clearly be- 
yond the scope of the present approach to account for the 
complete phase diagram of real manganites, which con- 
tains a great variety of charge, spin and orbital ordered 
states. The focus is on the doping region 0.15 < x < 0.5, 
where at low temperatures a band-description seems to 
be appropriate. In this regime, the phase diagrams of 
Fig. 4, determined for two characteristic parameter sets, 
(£1 = -0.125 eV, E 2 = -0.25 eV, t = 0.3 eV) and 
{Ei = -0.15 eV, E 2 = -0.175 eV, t = 0.3 eV), describe 
the major features of the metal-insulator transition lines 
in Lai_ x Ca x Mn03 and Lai_ x Sr x MnC>3 |35|. Further- 
more, even the absolute values of the critical concentra- 
tion x c (T = 0) ~ 0.17 and of the transition temperatures 
T c agree surprisingly well with the experimental data. 

As can be seen from the inset, the low-temperature 
phase contains besides the ferromagnetic metallic region 
finite domains of the localized polaronic phase (p^ < 1, 
and accordingly p^ > 0), the fraction of which increases 
with increasing temperature. Finally, at T c , there is an 
abrupt spill-over of holes between the delocalized and 
localized phases, which drives the metal insulator transi- 
tion. Since the conductivity is proportional to the num- 

ber of delocalized holes N¥ , a dramatic change at T c 
results. In the insulating high-temperature phase, the 
band description breaks down, and the transport proper- 
ties are dominated by incoherent small polaron hopping 
processes. For doping levels x < x c , the system behaves 
as a short-range correlated Jahn- Teller insulator. 
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IV. SUMMARY 

Although the theoretical understanding of the 
CMR phenomenon is still incomplete, double-exchange, 
electron-phonon and orbital effects are commonly ac- 
cepted as the main ingredients. Based on an in- 
depth analysis of important new experimental informa- 
tion about the inhomogeneous microscopic structure of 
the ferromagnetic metallic state, we come to the conclu- 
sion that in addition mixed-phase tendencies and percola- 
tive behavior play an important role in doped mangan- 
ites, in particular in the vicinity of the metal-insulator 
transition. To simulate these effects, we have studied 
a simple semi-phenomenological model, describing co- 
existing polaronic insulating and double-exchange domi- 
nated metallic phases with equal hole density. The frac- 
tion of localized and delocalized states was determined 
sclf-consistently. Below the transition temperatue T c , 
we found polaronic inclusions embedded in a dominant 
macroscopic metallic phase. Our approach, which is dis- 
tinct from that of the electronic phase separation concept 
involving states with different carrier density, provides an 
alternative explanation of the ferromagnetic metal to po- 
laronic insulator transition. The abrupt change, revealed 
in various electrical and magnetic properties at T c , was 
attributed to a collapse of the Zener state mainly caused 
by a percolative feedback mechanism. At T = the tran- 
sition is driven by doping and occurs at x c ~ 0.15 — 0.18. 
At finite temperatures, disorder due to intrinsic inhomo- 
geneities and magnetic scattering act in combination to 
reduce the mobility of the charge carriers. The calculated 
values of T c agree fairly well with the experimental ones. 
In conclusion, we believe that a further elaboration of the 
proposed percolative mixed-phase scenario, e.g., towards 
a direct calculation of the transport properties, may help 
to clarify many puzzling aspects of the CMR compounds. 
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APPENDIX: PERCOLATIVE PICTURE 

To support the assumption that the bandwidth of the 
Zener state depends approximately linear on the fraction 
of the ferromagnetic region (see Sec. II B, Eq. (j|)), let us 
consider a site percolation model on a finite hypercubic 
lattice with 64 3 sites and periodic boundary conditions. 
The lattice points are occupied with probability p. Adja- 
cent occupied sites will be connected by a hopping matrix 



element. The density of states of the resulting random 
tight-binding model, 



ft P = I>(cH+ c ] c l) 



(19) 



t p G {0, 1}, is determined numerically in two different 
ways, using kernel polynomial and maximum entropy 
methods On the one hand, we start from a single 

realization H. p and determine 400 Chebyshev moments of 
Hp using 100 random start vectors. On the other hand, 
we use 10 realizations of TL P and compute for each real- 
ization 200 moments from 50 random start vectors. The 
results basically agree. 
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FIG. 5. Density of states (DOS) (upper panel) and in- 
tegrated DOS (N(E)) for the tight-binding site percolation 
model with different probabilities p = 0.1,0.3,... ,1. The 
inset gives the bandwidth as a function of p. 

Figure | displays the (averaged) density of states for 
various values of p. Self-evidently, p now is simply a 
model parameter, i.e., there is no feedback effect as con- 
sidered in Sec. II B. A significant change of the spec- 
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trum around the site percolation threshold, p c ~ 0.31, 
is observed. Below p c , the DOS is mainly confined to 
the interval [-2,2], like for a one-dimensional model, and 
mostly consists of isolated peaks with different spec- 
tral weight. This can be seen more clearly from the 
stcplike behaviour of the integrated density of states, 
N(E) = J E w/2 DOS(E')dE', shown in the lower panel. 
The steps can be attributed to metallic "islands" (finite 
clusters) containing delocalized electrons. They form a 
cluster comparable to the system size first at p = p c . 
Above p c , N(E) becomes a continuous function and the 
DOS is evocative of that of a simple cubic tight-binding 
model but with a reduced bandwidth. That is, although 
the percolative cluster occupies most of the bulk, making 
the whole system metallic, smaller conducting and insu- 
lating regions are embedded into it. These enclaves act 
as small scattering centers, causing the band narrowing. 
The inset demonstrates the almost perfect linear depen- 
dence of W on p, in particular at large p, which in reality 
is the region of interest (cf. the inset of Fig. El). 
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